Spin dynamics in electron-doped pnictide superconductors 
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The doping dependence of spin excitations in Ba(Fei_iCo2,)2 As2 is studied based on a two- 
orbital model under RPA approximation. The interplay between the spin-density-wave (SDW) and 
superconductivity (SC) is considered in our calculation. Our results for the spin susceptibility are in 
good agreement with neutron scattering (NS) experiments in various doping ranges at temperatures 
(T) above and below the superconducting transition temperature Tc. For the overdoped sample 
where one of the two hole pockets around T point disappears according to ARPES, we show that the 
imaginary part of the spin susceptibility in both SC and normal phases exhibits a gap-like behavior. 
This feature is consistent with the "pseudogap" as observed by recent NMR and NS experiments. 

PACS numbers: 74.70.Xa, 74.25.Ha, 75.30.Ds 
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The recent discovery of the iron arsenide supercon- 
ductors whose parent compounds exhibit long-range 
antiferromagnetic (AF) or spin-density-wave (SDW) or- 
der similar to the cuprates provides another promis- 
ing group of materials for studying the interplay be- 
tween magnetism and superconductivity (SC). Espe- 
cially, the electron-dj3ped pnictide superconductors like 
Ba(Fei_a;Co2;)2As2 [Si] has emerged as one of the most 
important systems due to the availability of large ho- 
mogeneous single crystals. The phase diagram [4|-|6| for 
these materials indicates that the parent compound upon 
cooling through Tat ~ 140K develops a static SDW 
order. Increasing the doping of Co, the SDW order is 
suppressed and the SC order emerges as the temperature 
(T) falls below T^. The SDW and SC orders coexist in 
the underdoped samples By further increasing 

the Co concentration to the optimally doped regime, the 
SDW order disappears. These experimental results pro- 
vide compelling evidences for strong competition between 
the SDW and SC orders. 

Recently, several neutron scattering (NS) experiments 
have been carried out to probe the spin dynamics in these 
materials [1, @, Isl-flil , and the spin excitation spectrum 
was fitted by usiiig a Ji , J2 Heisenberg model based on 
localized spins However, while the parent com- 

pounds of the cuprates are Mott insulators with large 
in-plane exchange interactions [l6| . the parent iron ar- 
senides are bad metals and remain itinerant at all dop- 
ing levels. Magnetism in these materials are most likely 
to originate from itinerant electrons and the AF order 
is a result of SDW instability due to Fermi-surface nest- 
ing Theoretically, at present, the variation of the 
spin susceptibility with doping remains less explored. 
The spin susceptibilities were mostly studied in the op- 
timally doped compounds without SDW [3| or in the 
parent compound without SC [19]. In this work, we 
adopt Fermi-liquid mean field (MF) theory to study the 
static SDW and SC, and employ the random-phase ap- 



proximation (RPA) to investigate the spin dynamics in 
Ba(Fei_a;Coa:)2As2 from the imaginary part of the dy- 
namic spin susceptibility. It is expected that the present 
approach is much more justified than the localized model 
for examining the spin fluctuations in the iron arsenides, 
as suggested by both experiments [T^, and theories 
[isl [l9| . We show that the calculated spin susceptibil- 
ities are in qualitative agreement with several NS and 
NMR experiments in various doping ranges. 

We start with a two-orbital model by taking into ac- 
count two Fe ions per unit ce ll 1 2011 . The reason we 
adopt this model is its ability [21| to qualitatively ac- 
count for the doping evolution of the Fermi surface 
and the asymmetry in the SC coherent peaks as ob- 
served by the angle resolved photo-emission spectroscopy 
(ARPE S) |22l| and the scanning tunneling microscopy 
(STM) [23| experiments on Ba(Fei_j;Coa;)2As2. The 
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Hamiltonian of our system can be expressed as 
H = Ho + Ha + Hint- Hq is the tight-binding Hamil- 
tonian and can be written as Hq = J2wa V'kCT-^'^kV'ka 

.... 



2M, where V'L = (cAo.k<T' c^nw^, 
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the creation operator with spin a = —1 (J,) or 1 (t), in 
the orbitals (0,1) = {dxz,dyz) at the sublattice A {B), 
and 
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where eam = ~2(t2COskx + 1 3 cos ky 
—2{t2 cos ky + t^coskx), Sxy.yi = —2t4{coskx + cosk^), 
and ET.k = —4^1 cos ^ cos ii_4 are the hopping pa- 
rameters and n is the chemical potential. Throughout 
the paper, the momentum k is deflned in the tetragonal 
notation. 



The pairing term is Ha = Eksal^kC^ 



sa,kt sa,— k^ 



h.c). Here, we assume there exits only next-nearest- 
neighbor intraorbital pairing with extended s— wave pair- 
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ing symmetry, and the SC order parameter is Ak 



0.018 
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2° (cos kf r ~\~ COS kii) similar to that obtained by spin fluc- 



tuations 



17[, where 



Ao = — y^(cosA;2; + cosfcy)(c 
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with Vnnn being the attractive pairing interaction. 

Hint is the on-site interaction term which includes the 
Coulombic interaction and Hund coupling Jh, following 
Refs. 3 21, 24 [, it can be expressed as 



Hint = U ''^nsa,ij^nsa,iji + {U 



ijs 

— 2 J^f Sso.ij • Ssi^ij, (3) 



where {i,j} denotes the unit cell, s = (A) or 1 (B) is 
the sublattice index, and a = (dxz) or 1 (dyz) represents 
the orbital, risa.ija and Ssa.ij are the density and spin 
operators in the orbital a at the sublattice s of the unit 
cell respectively. According to symmetry, we have 

U' = U - 2Jh [11] ■ In the MF approach, we linearize 
Hint in momentum space as (21I . [25| 

- y(C^ + ^H) '^^L,k+Q^C3„,ko-, (4) 

ksao" 

where n = 2 -I- x is the number of electrons per lattice 
site and Q = (7r,7r). The SDW order parameter is to = 
Lk<T ^("^L.k+QtrCsQ.kff), with N being the number of 
unit cells. 

The effective MF Hamiltonian is then given by 
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where (^j, = (V'^^, i/'k+qt, V'-ki, V'-(k+Q);), A/k = ^/k + 
f (3i7 - 5Ji/)/, and R = -f{U + Jh)L / is a 4 x 4 

unit matrix and ^j^. means the summation extends over 
the magnetic Brillouin zone (MBZ): — tt < kx i ky ^ 
TT. The MF Green's function matrix can be written as 
.g(k,T) = -{TrMr)^liO)) and g{KiPn) = ^kW^k^, 
where W^kij ~ ^^i ^ ^ki ) ^ ^ and ^k is a unitary matrix 

that satisfies (A]{^Wk^k)ij = (^ijAki. 

First, we solve the MF equations self-consistently to 
obtain m, Aq and fi at different doping levels x and 
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FIG. 1: (Color online) (a) m? as a function of T close to Tc 
at a; = 0.06. Inset shows the calculated phase diagram, (b) 
Imx^^^ici,^) at q = (tTjTt) as a function of oj by changing 
iijJn to OJ + irj, at X = and different temperatures T. (c) 
Imx^^^{<i,oj) as a function of energy transfer uj and mo- 
mentum q, at a:: = and T = 0. The momentum is scanned 
along (qx/TTjqy/Tv) — {h,h). (d) Constant-energy scans along 
the {qx/iT, Qy/i^) ~ (h, h) direction at a; = and T = Q. Suc- 
cessive cuts are displaced vertically for clarity. The energy 
transfer is a; = 0.09 (black solid), 0.4 (red dash), 0.93 (green 
dot), 1.2 (blue dash dot), and 2 (cyan short dash). The damp- 
ing rate r) — 0.04. 



temperatures T. The magnitudes of the parameters 
are chosen as ti_4 = 1,0.4,-2,0.04 U = 3.4, 

Jh = 1-3, Vnnn — —1-2 f2l'| , and the number of unit 
cells is 257 x 257. Throughout the paper, the energies 
are measured in units of \ti\. The calculated phase di- 
agram as shown in the inset of Fig. HJa) reproduces 
the result based on Bogoliubov-de Gennes (BdG) equa- 
tions (2l| and is also consistent with the experiments on 
Ba(Fei_^Co:r)2As2 Here the SDW and SC are 

competing with each other. If there is no SDW, SC 
would show up even in the parent compound. The pres- 
ence of SC also suppresses SDW. For example, in the 
underdoped [x = 0.06) compound with Tjv ~ 0.16 and 
Tc ~ 0.05, the calculated magnitude of (proportional 
to the magnetic Bragg peak intensity) at T = is re- 
duced by 4% relative to that of the maximum intensity 
at Tc [see Fig. [Ija)], and this result is consistent with the 
neutron diffraction experiments 0, Q • 

Then, we investigate the spin dynamics in 
Ba(Fei_^C0:r)2As2 for a: = 0, 0.06, 0.1, and 0.2, 
corresponding to the undoped, underdoped, opti- 
mally doped and overdoped compounds, respectively. 
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FIG. 2: (Color online) (a) /mx^^'^(q, t^) at q = (tt,tt) as a 
function of uj, ai x = 0.06 and different temperatures T. (b) 
/mx^^'*(q, as a function of energy transfer ui and momen- 
tum q, at X = 0.06 and T — 0. The momentum is scanned 
along {qx / , Qy / tt) = {h,h). (c) and (d) ((e) and (f)) are sim- 
ilar to (a) and (b), respectively, but at a; = 0.1 (x = 0.2). 
Inset in (e) shows the x = 0.3 case. 



The MF spin susceptibility is 

s 



ra,sf3{0)^ ' ■ x 



^ J2[°j)=(^% [Pim.nj {q)+Pi+4m,nj+4:iq)] , here, r, s, t, u 
label the sublattice indices, a, /3, 7, (5 represent the or- 
bitals, TO = 2r + a + l,n = 2s + /3+ 1,0 = 2t + 7+ l,p = 
2u + 5 + l, and Pim,nj{q) ^ ~Jn I^k,j9„ 9im{k)g„jik + q). 
Here we used k = (k, ipn) and g = (q, iw„). 

We then use RPA to take into account the residual 
fluctuation of Hint beyond MF. The RPA spin suscepti- 
bility is determined by the matrix equation x^^'^io) — 
E.to^{x^'?)[/-rx°(g)]-H™, where / is a 16x 16 unit 
matrix and the nonzero elements of the interaction vertex 
are: for a = /3 = 7 = ^, r™'^f = 2t7; for a = /3 7^ 7 = J 

or a = 7 ^ /5 = <5, r™_vt = 2 J^. 

In the parent (x = 0) compound, the RPA spin sus- 
ceptibility [see Fig. W[h)] in the paramagnetic state at 
T ~ 0.27 (Tjv ~ 0.25) shows a linear energy depen- 
dence for a; < 0.1, suggesting gapless excitations lol 14|. 
On the other hand, in the SDW state at T = 0, the 



spin excitation intensity is close to zero below uj w 0.13, 
similar to a spin gap 0, However, the gap is not 
sharp since a sharp gap would produce a stepwise in- 
crease in intensity at the gap energy which is unlike the 
more gradual increase seen here. Figures [TJc) and [TJd) 
show Imx^^^{<l,oj) as a function of energy transfer oj 
and momentum q along (fix I , % I t^) = {h,h) direction 
at T = 0. As we can see, there is almost no detectable 
intensity below uj w 0.13, again illustrating the opening 
of the spin gap. The excitations are peaked at Q = (tt, tt) 
and UJ w 0.93, at higher energies, the response is seen to 
split and broaden due to the dispersion of the spin waves. 
By tracking the peak positions in Fig. [UJc) , the spin- wave 
disp ersion relation can be fitted as ujq = y^A^ + iPq^ 
fwL 15, 2^, where A « 0.93 is an energy gap, v « 9.02 
is the spin-wave velocity, and q is the reduced wave vec- 
tor away from (1, 1) along the (h, h) direction, consistent 
with the experimental observations [l3|- The origin 
of the spin gap can be understood in terms of the Fermi 
surface at a; = as shown in Fig. 2(a) in Ref. 2l|. In 
the paramagnetic state, large parts of the two hole pock- 
ets around F = (0, 0) and two electron pockets around 
M = (7r,7r) are nested by momentum (7r,7r), thus giv- 
ing rise to the gapless excitations at T = 0.27. But at 
T — 0, the SDW order will gap most parts of the original 
Fermi surface, leaving only tiny ungapped Fermi surface 
pockets connected by (tt, tt) along the F — M line, so in 
this case, for small energies, the imaginary part of the 
MF spin susceptibility is close to zero, while its real part 
does not fulfill the resonance condition, leading to a spin 
gap opening in the RPA spin susceptibility. 

The RPA spin susceptibility [Fig. dja)] at a; = 0.06 
suggests that the excitations above Tc are gapless, al- 
though the intensity is very small at low energy. This 
may be the reason why above Ref. [1] claims the 
excitations are gapless while Ref. [§] concludes they 
are gapped. Below Tc, the intensities below uj w 0.06 
and above uj « 0.36 are suppressed and the weight is 
transferred to form a resonance at w^es ~ 0.14. Since 
m{T = 0.09) « 0.129 and m{T = 0) « 0.131, our re- 
sults seem to agree with Ref. |^], which claims the res- 
onance is produced by suppressing low energy spectral 
weight, rather than Ref. [8|], where the spectral weight is 
considered to be transferred from the ordered magnetic 
moments. In addition, figure[2]Jb) shows that below ujres, 
commensurate spin excitation prevails, in agreement with 
experimental observation [8j], and it becomes incommen- 
surate when the energy is above ujj-esi notably between 
UJ « 0.15 and uj ~ 0.18, which we predict to be measur- 
able by NS experiment. The spin excitations can extend 
beyond uj = 1, with smeared out and broadened features 
for UJ > 0.2. 

At X = 0.1, the SDW order is completely suppressed, 
and SC emerges for T < Tc ~ 0.06. The excitation spec- 
trum [Fig. m^c)] shows that in the superconducting state 
at T = 0, a gap below uj ss 0.08 develops and there is a 
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resonance above the gap energy peaking at ujres ~ 0.24, 
in agreement with the NS ejcperiments on the optimally 
doped Ba(Fei„2;Co2;)2As2 _11H14|. Furthermore, figure 
^d) shows that the spin excitation is incommensurate at 
low energy (0.18 < w < ojres) which still need to be veri- 
fied by experiments, then it switches to a commensurate 
behavior between uj ~ ujres and lo w 0.3, and becomes 
broad at higher energy, consistent with Refs. 12. 13. 15|. 
In the normal state at T = 0.12, the spectrum is replaced 
by broad gapless excitations with a linear energy depen- 
dence for w < 0.1 [11,1 14| . We notice a marked similarity 
between the spin excitations in the normal state of the 
optimally doped compound and those in the paramag- 
netic state of the parent compound as observed in Ref. 
[3], suggesting a common origin of spin fluctuations in 
both of them. 

In contrast, the spin excitations in the overdoped 
{x ~ 0.2) compound [Fig. EJe)] show gap-like behav- 
ior in both the normal and superconducting states. The 
origin of the gap may be due to one of the two hole pock- 
ets around F vanishes and the other one shrinks dramati- 
cally in the overdoped re gion according to ARPES exper- 
iments 22 1 and theories 2^, 21|. Under such a case, due 
to the lack of interband scattering between the hole and 
electron pockets, the imaginary part of the spin suscep- 
tibility is strongl y su ppressed and gives rise to the pseu- 
dog ap behavior 27| which has been observed in NMR 
[28j and NS [13] experiments in the electron overdoped 
Ba(Fei_a;Coa;)2As2. But in Ref. [27|, the pseudogap is 
associated with the vanishing of one of three hole pock- 
ets around F, where experimentally there is only one 
hole pocket at this doping level as observed by ARPES 
[2^ . The spin excitations in the superconducting state 
at r = [Fig. ^i)] are broader and weaker than those 
in the underdoped and optimally doped compounds, sug- 
gesting the importance of the hole pocket in enhancing 
the spin fluctuations. 

At X = 0.3, both the two hole pockets around F disap- 
pear 2^, 21 1, our calculations show that SC is completely 
suppressed and the spin fluctuations are extremely small 
[The inset of Fig. [UJe)]. This further indicates the cor- 
relation between the electronic band structure and mag- 
netism, and supports the scenario that the spin fluctua- 
tions in the underdoped regime, which serve as a precur- 
sor to SC, originate from quasiparticle scattering across 
the electron and hole pockets. 

In summary, we have systematically investi- 
gated the doping dependence of spin excitations in 
Ba(Fei_a;Cox)2As2, ranging from the parent to over- 
doped regime. In the parent compound, the spin 
excitations are gapless in the paramagnetic state and 
become strongly suppressed at low energy in the SDW 
state due to the opening of gaps on most parts of the 
original Fermi surface. For underdoped and optimally 
doped samples, the spin gaps and resonances at (tt, tt) 
only occur in the SC state. On the other hand, the spin 



excitations in the overdoped compound show gap-like 
behavior in both the normal and SC states due to the 
vanishing of one hole pocket around F, leading to a 
"pseudogap" behavior at this doping level. All the 
obtained results are in qualitative agreement with exper- 
iments. The changes in the spin dynamics at different 
doping levels may reflect changes in the electronic band 
structure and suggest a strong correlation between SC 
and magnetism. 
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